Preprint typeset in JHEP style - HYPER VERSION 



EFI 10-32 



Low-energy Observables and General Gauge 
Mediation in the MSSM and NMSSM. 



Arun M. Thalapillil 

Enrico Fermi Institute and Department of Physics, 
University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637 
E-mail: madhavSuchicago . edu 

Abstract: We study constraints on the general gauge mediation (GGM) parameter space 
arising from low-energy observables in the MSSM and NMSSM. Specifically, we look at 
the dependence of the spectra and observables on the correlation function ratios in the 
hidden sector where supersymmetry is presumably broken. Since these ratios are not a 
priori constrained by theory, current results from the muon anomalous magnetic moment 
and flavor physics can potentially provide valuable intuition about allowed possibilities. It 
is found that the muon anomalous magnetic moment and flavor-physics observables place 
significant constraints on the GGM parameter space with distinct dependences on the 
hidden sector correlation function ratios. The particle spectra arising in GGM, with the 
possibility of different correlation function ratios, is contrasted with common intuition from 
regular gauge mediation (RGM) schemes (where the ratios are always fixed) . Comments are 
made on precision gauge coupling unification, topography of the NLSP space, correlations 
of the muon anomalous magnetic moment with other observables, and approximate scaling 
relations in sparticle masses with respect to the high-scale correlation function ratios. 
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1. Introduction 

Supersymmetry (SUSY) is considered to be one of the most promising extensions to the 
standard model (SM), since it potentially explains many of the problems in it. In all 
supersymmetric extensions of the SM it is required, for reasons of viability, that SUSY 
be broken in a hidden sector which then gets communicated to the visible sector (see for 
example [1] and references therein). An appealing implementation of this SUSY breaking 
paradigm has been regular gauge mediation (RGM) [2] since it solves the flavor problem 
and at the same time is parsimonious compared to, for instance, gravity mediation. The 
phenomenology of the RGM implementations have been studied extensively in the MSSM 
and its extensions (see for example [2] and related citations). 

Recently it was pointed out in [3] that the parameter space of gauge mediation is 
potentially larger and a convenient parametrization to account for and study this was 
introduced there with the terminology of general gauge mediation (GGM). Since then 
there have been many studies exploring the phenomenology of this enlarged parameter 
space [4, 5, 6, 7, 8, 9, 10]. Also, on the theoretical side, the ideas have since been extended 
[11, 12, 13, 14, 15, 16]. The mass spectra in GGM was first investigated in [4] where it was 
pointed out that the spectra could be very distinct from minimal gauge mediation. For 
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the MSSM, an investigation of the NLSP topography for various GGM correlation function 
ratios was instigated in [5] and more recently the prompt decays of NLSPs and co-NLSPs 
were investigated in [6] and [7] . Low-energy observables with fixed ratios of the correlation 
functions were investigated in [8] and in [9] benchmarks points in GGM were discussed 
(again for fixed ratios), in the context of LHC searches. 

We wish to complement the studies of [4, 5, 8, 9], by exploring the effects of dif- 
ferent hidden sector correlation function ratios on low-energy observables and the NLSP 
topography, in the context of the Minimal Supersymmetric Standard Model (MSSM) and 
the Next-to- Minimal Supersymmetric Standard Model (NMSSM). The hope in the present 
study is that since these ratios are not a priori constrained by theory, current results from 
the muon anomalous magnetic moment and flavor-physics observables can be useful guides. 
If low-energy SUSY is discovered at the LHC and if the soft masses can be measured to very 
good accuracy (~ 1% uncertainty) sometime in the future, then it may even be possible to 
determine the correlation functions at the high scale to good precision using renormaliza- 
tion group invariants [17, 18]. More recently, the role of the GGM messenger scale in the 
context of mass sum rules and RG invariants was considered in [19, 20] 

In the next two subsections we introduce our notations and definitions for the MSSM 
and NMSSM, and then proceed to review relevant features of RGM and GGM that are the 
main focus of our investigation. Then in section 2 we briefly review some of the low-energy 
observables and collider bounds that we use to constrain GGM and to explore correlation 
function dependences. In section 3 we present our numerical study along with results and 
observations. Section 4 is the summary. 

1.1 The MSSM and NMSSM 

MSSM is a supersymmetric extension of the SM where the effects of SUSY breaking are 
parametrized by super- renormalizable soft terms. For reasons of anomaly cancellation and 
holomorphicity of the superpotential there are two Higgs doublets in the MSSM apart from 
the sparticles. In this subsection we fix our notation. 

The superfields (denoted by a ' " ' ) in the MSSM are defined as 



where all the generation indices have been suppressed and the respective SU (3)c x SU (2) l x 
U{1)y assignments are shown in brackets. The MSSM superpotential (with all the gener- 
ation indices again suppressed) is given by 




(1.1) 



MSSM — 



YuHu-QU'R + YdHd-Q D'ji + Y.Hd-LE^ + fiHu-Hd. 



(1.2) 
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where we define A - B = eijA^BK The corresponding soft SUSY breaking terms in the 
MSSM may be parametrized as 

-A^r = BB + M2WW + M3 ~g~9 + h.c.) + mljHu\^ + mjjjHd]^ 

+m||Q2| + ml\U%\ + m|,|D|| + m| jL^] + m||£;|| 
+iYuAu Q-HuU'k- YdAd Q-HdDji- Y^Ae L-HdE'ji + h.c.) 
+{BiiHu-Hd + h.c.) , (1.3) 

where B,W and g are the gauginos corresponding to U{1)y,SU{2)l and SU{3)c respec- 
tively. We have the usual definition for tan (3 

tan/3 = — , (1.4) 

Vd 

and the soft mass terms, j4-terms and Bfi are defined in the standard way [1]. The MSSM 
Higgs potential has four free parameters 

H,Bfi,m'jj^,m'jj^ . (1.5) 

Using the minimization conditions we may trade some of these parameters for others, for 
instance tan /3 and M^o . 

Augmenting the MSSM with a singlet chiral superfield N we may define the NMSSM 
superpotential as (see for example [21]) 

Wkmssm = Wmssm + \NHu-Hd + (fN + ^/i'iV^ + ^iV3 . (1.6) 

The corresponding soft SUSY breaking masses and couplings in the NMSSM are 

_^NMssM ^ BB + M2WW + M3M + h.c.) + m%jHu\^ + mljHdl^ + m%\N\^ 

-fm^lQ^I + m||C/^| + m|,|D|| + m||L2| + m||£;|| 
+{YuAu Q-HuUfi- YdAd Q-HdD'ji- Y^A^ L • E'j^ + h.c.) 
+iXAx Hu-HdN+ ^kA, + mlHu-Hd + ]^m% + N + h.c.) .(1.7) 

Usually a simpler NMSSM potential is considered by imposing scale invariance whereby 
the parameters //, /i', .^i;', mg, and ^at are vanishing. This potential has a discrete Z3 
symmetry and the Higgs sector is described by seven parameters 

X,K,Ax,A^,mjj^,mH^,m%. (1.8) 

One may again use the minimization conditions for the NMSSM Higgs potential to trade 
some of these for other quantities. 
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1.2 Regular Gauge Mediation and General Gauge Mediation 

The need for a separate SUSY breaking sector to accommodate viable phenomenology is 
well known. In this paradigm the SUSY breaking is expected to occur in a separate hidden 
sector and the SUSY breaking is then communicated to the visible sector directly /indirectly. 
If the mediators between the hidden and visible sector are physical particles then they 
constitute the messengers of the SUSY breaking. The mass scale associated with these 
messenger superfields is termed the messenger scale. A particular realization of this scheme 
is gauge mediation where the messenger fields are assumed to be flavor blind. 

In regular gauge mediation (RGM), the gaugino and sfermion soft masses arise from 
loops involving messenger fields and are given by [2] 



F 



M„ 



where 



fix) = ^ [(1 + x) log(l + x) + {l-x) log(l - x)] , 

g{x) = ^ [log(l + x) - 2Li2(x/[l + x]) + iLi2(2x/[l + x])] + (x ^ -x) ; (1.10) 



X = F/M"^, C2(/|r) are the quadratic Casimirs, Mm is the messenger scale and Nm is the 
number of copies of the messenger particles in the loop (or in other words the sum of the 
Dynkin indices). The RGM includes both minimal gauge mediation, where Nm = 1, and 
non- minimal gauge mediation where Nm > 1- 

For small x, f{x) — >• 1 and g{x) — )• 1 in Eq. (1.9). This implies that for Nm = 1 

Mr ^ m-j , (1.11) 

since the sfermion mass squared is at 2-loop while the gaugino mass is at 1-loop. In RGM 
we can certainly tune the gaugino to sfermion mass ratio to some extent by increasing Nm, 
but beyond a point it is very difficult without further inputs to get a viable model along 
with gauge coupling unification. We will see that in GGM it is easier to get large or small 
gaugino to sfermion mass ratios naturally. 

In RGM, the gravitino is always the LSP since 

F 

misp = rn3/2 ~ — , (1.12) 

and Mm ^ Mpi. At low (electroweak) scale ai,a2 < as and hence the NLSP is usually 
the X or I in large regions of the parameter space. 

Moreover due to gauge coupling unification and the fact that Ma/ ga is an RG invariant 
(to one-loop order) it is always true, at any scale, that 

Mi_M2_M3 

91 9i ai 
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at this order. 

The corresponding expressions in general gauge mediation (GGM) are [3] 

Mr = g^rMs0\o) , 

= gfYfC + E gtC2{f\r)M^Ar , 

J r=l 

where 



(1.14) 



(1.15) 



Sy^(O), cP (where pis a. vector, fermion or scalar index) are associated with the current- 
current correlators in the hidden sector, C2(/|r) are again the quadratic Casimirs, is a 
possible Fayet-Iliopoulos term (D-term) and Mg is some characteristic soft SUSY-breaking 
scale associated with the hidden sector. With this definition note that the correlation 

functions {By^Ar) are dimensionless. Also note that the subscript 'r' labels the associated 

~ 1/2 

SM gauge group. Henceforth we will abbreviate Br (0) as just Br, not to be confused 
with the notation for a bino [B) that has no subscript index. 

In the GGM parametrization the interpretation of the messenger sector is enlarged 
to accommodate more general scenarios, for instance, cases where the sector is strongly 
coupled and there are no explicit messenger particles. Note that the RGM limit can be 
obtained from the GGM scenario by considering points in the vicinity of 




Br = \h^, (1.16) 



for the case when there is no inter or intra hierarchies among the Br and Ar. By considering 
points near the above region we can accommodate both minimal and non-minimal RGM 
cases. 

Thus now we have the varied possibilities 

Mr w mj, Mr <. nij, Mr > mj , (1-17) 

depending on the correlation function hierarchies (i.e. hierarchies among Br and Ar). Also 
note that it is the case that in general now 

9i 92 9i 

at any scale, even when there is gauge coupling unification, due to possibly different Br = 

On a final note, the Fayet-Iliopoulos term [Q that makes an appearence in the ex- 
pression for the sfermion masses is potentially dangerous since it could lead to tachyonic 
masses. It is usually set to zero due to this by imposing some discrete symmetry like 
messenger parity [3]. 



- 5 - 



2. Low-energy observables and constraints 
2.1 Anomalous magnetic moment of the muon 

The anomalous magnetic moment of the muon has been measured to very good precision. 



Current measurements yield a value [23] 

[a^]^^^ = (11 659 2080.0 ± 54stat ± 33sy.t) x 10 
The theoretical contributions may be divided as [22, 23, 24, 25] 
[a,,]^^^ = (11 658 4718.1 ± 1.6) x 10^^^ , 
a 



11 



(2.1) 



a 



-"MJHad. LO ~ 
M-lnad. NLO 



(6955 ± 40exp ± 7qcd) x 10^ 
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(-97.9 ± O.Scxp ± 0.3rad) X 10- 



-11 



;i05±26) X 10^ 



Jcxp 

-11 



(2.2) 



(154 ± lhad ± 2Higgs) X 10- 



-11 



-11 



Klxot. SM = (11 659 1834 ± 4lLOhad ± 26NLOhad ± 2othcr) X 10- . 
In the above 'had.' means hadronic and 'LBL' stands for light-by- light scattering. For 
the hadronic vacuum polarization corrections we have quoted the value from a{e~^e~ — )• 
Hadrons) low-energy data. Comparing the theoretical prediction and the most current 
experimental result gives a discrepancy 



[Aa^ 



Exp. 



(246 ± 80) X 10 



-11 



(2.3) 





Figure 1: The well known 1-loop contributions to (17—2)^1 in the MSSM. The dominant contribution 
comes from the chargino diagram. For brevity the diagrams are shown in terms of the mass 
eigenstates. 

In the MSSM the 1-loop contributions to the muon anomalous magnetic moment come 
from loop diagrams with a chargino-sneutrino and neutralino-smuon in the loop and are 
given by [26] 



,SUSY,1L 



X 



(2.4) 
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with 



% = 


167r2 ^ 

k . 


x° 


167r2 ^ 

i,m _ 



2^77, -]- 



12m 



3mi 



l-ln 



(2.5) 
,(2.6) 



where i = 1 . . . 4 and k = 1,2 denote the xo ^^nd x-t, m = 1, 2 denotes the jl index. Also 
we have the definitions 



Qk 



Pim 



Pirn 



-52Vfci , 

yijUk2 , 

V2g,NaU^^ + y^N,sU^,. 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



The matrices U, V, N, relate the mass eigenstates to the interaction eigenstates for the 
gauginos, Higgsinos and sfermions and are defined as 



Xa = Uabi'h ' 
Xa = Nab^P^b , 
fa = Ulf^fh ■ 



(2.11) 



where the LHS states are mass eigenstates (charginos, neutralinos and sfermions) and h 
labels the chirality state. With the definitions Xim = ''^'^o/'mfi^ and Xk = m^^/mf,^ the 
loop functions are [26] 



Xiix) = 


(1 


2 

— X 


^ [2 + 3x — 6x^ + x^ + 6x log x] , 


X2{x) = 


(1 


3 

— X 


g- [ — 3 + 4x — — 2 log x] , 


Yi{x) = 


(1 


2 

— X 


^ [l — 6x + 3x^ + 2x^ — 6x^ log X 


Y2{x) = 


(1 


3 

— X 


-g [l — + 2x log x] . 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



To get some intuition about the above 1-loop contributions let us make some simpli- 
fying assumptions. Let us take tan/3 to be large with all soft masses roughly equal to a 
common scale MgusY] then the above expressions may be approximated as 



X" 9l - 9l 



/7A. 

" - 1927r2 M? 



sign(/iAf2) tan/3 , 



7X 



9l 



SUSY 
2 



m 



(2.16) 
(2.17) 
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where real parameters and equal signs of Mi and M2 have been assumed. From the above 
approximate expressions it is clear that the contribution to is numerically larger than 
the X contribution. It is also noted that the sign of the SUSY contribution is determined 
by the sign of /i. 

We can estimate a large class of 2-loop contributions by noting that the large logarithms 
arising from QED corrections to one-loop SUSY diagrams may be quantified as 

^^SUSY+QED ^ ^^l-Loop{SUSY) ^SUSY 



(2.18) 



This leads to a reduction by a few percent of the LO contributions. 

The contributions to at 1-loop in NMSSM are identical to those in the MSSM 
provided that we take into account the additional singlino state in the x calculation. The 
CP-even and CP-odd parts of the singlet field in the NMSSM will mix with the neutral 
components of the Higgs doublets. In certain regions of the parameter space it is therefore 
possible to have a relatively light CP-odd Higgs scalar which can provide a significant 
contribution to a^. 

In the SM or MSSM, Higgs effects are usually negligibly small because of the current 
lower bounds on the Higgs masses. The SM 1-loop Higgs diagram, for instance, is about 
2-3 orders of magnitude below experimental sensitivity for nih ^114 GeV. 

In the NMSSM on the other hand the mass bounds are relaxed relative to the MSSM. 
In the NMSSM the lightest CP-odd Higgs boson ai can be as light as a few GeV and still 
satisfy constraints from flavor physics, especially from Br{Bs — )• for low values of 

tan/3 or when the loop-induced coupling is small. 

The SUSY Higgs contributions (CP-even, CP-odd and charged) to may be calcu- 
lated as [27] 

e^lL CP-even ^f^'^fi ST^ ^2 ^^(^ ~ ^) dx 

4 V 27r2 cos2 /3 ^ Jo 2^('^'Hy n 



= - ""'"^"'^ V /' , (2.20) 

"^-^ A . - 1 + i^^) 

where Xij^Yij are the Higgs mixing matrices for the CP-even and CP-odd cases. 

It has been noticed previously that the 1-loop and 2-loop Higgs contributions come 
with opposite signs. Moreover it is seen from the above formula that the 1-loop CP-odd 
scalar contribution is negative. This will be potentially relevant to the NMSSM case where 
there is an extra CP-odd Higgs that could be relatively light. 

2.2 Flavor-physics constraints 

Flavor physics can place strong constraints on supersymmetry since flavor-changing (FC) 
contributions that are in general loop-suppressed can become important due to tan /3 en- 
hancement. We use the flavor observables shown in Table 1 in our study of the COM 
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Flavor-physics Observables 


IVIeasurements 


2(7 bounds 


iiT[n — > As7j 


(o.O^ ih U.zo ih U.Uyj X iU 


z.io X iu s Jjr(^o — > ) s 4.oy x lu 


DV[ns fJ- fJ- ) 


< o.o X iU 


dT[JDs — f H jl j < D.D X iU 


A ]\/r 


{ ^ -\- Ci c\A\ \y in^l 
(O.U( ±U.U4j X iU 


4.yy X iu < L\ivi(i < o.io x lu 




(1.777 ±0.012) X 10+1 


1.753 X 10+1 < AA/, < 1.801 x 10+^ 


Ao{B ^ K*7) 


(3.1 ±2.3) X 10-2 


-1.7 X 10-2 < Ao < 8.9 X 10-2 


Br(i^ ^ /izy) 

r^ylTT V 111/ I 


0.6358 ±0.0011 
1.004 ±0.007 


0.6257 < ^ , i < 0.6459 
0.990 < i?f23 < 1-018 




(1.41 ±0.43) X 10"'' 
1.28 ±0.38 


0.39 X IQ-'^ < Bt{Bu TVr) < 2.42 x 10""^ 
0.52 < Rri,^ < 2.04 




(8. 6 ±2. 4 ±1.1 ±0.6) X 10~3 
0.416 ±0.117 ±0.052 


2.9 X 10-3 < Bi{B D^TVr) < 14.2 x 10-^ 
0.151 < ^De^ < 0.681 


Br(Ds — > TVr) 

Br(L», ^ /izv^) 


(5.7 ±0.4) X 10-2 
5.8 ±0.4 X 10-3 


4.8 X 10-2 < Bt{Ds TVr) < 6.6 X 10-2 

4.9 X 10-3 < Bi{Ds fiv^) < 6.7 x 10-^ 



Table 1: Flavor-physics observables [29, 30, 31, 32, 33, 34, 35]. Most of the intervals are adopted 
from the relevant tabic in the Superlso manual [36]. For details on the derivation of these intervals 
the reader is pointed to references therein. The flavor observables listed in the table are defined in 
the text below. 



parameter space. The idea is that the hidden sector correlation functions {Br, Ar) deter- 
mine the sparticle spectra at some large scale which then RG evolve down to the low scale 
and potentially contribute to the Wilson coefficients of the various meson transitions. 

The rare branching ratios Br{B — ?• -^^57), Br{Bs — > fi'^fi"), Br{B — t- tv) and Br{B — )• 
Dtu) are used with their 2cr intervals to check for viability of points. Note that the best 
bounds on a charged Higgs in the MSSM are from Br[B — )• ^^7) as opposed to collider 
studies. This is a telling feature of the complementary nature of collider and low-energy 
constraints. 
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The b ^ s'j transitions may be parametrized by 

8 



n 



eff 



^ p=u,c i=l 



(2.22) 



where Oi{fi) are the relevant operators and Ci{^) are the Wilson coefficients evaluated 
at the scale ^ corresponding to these operators. They may be expressed in the standard 
operator basis [37] by: 



O2 
O3 

O7 



{s-f^T''pLc){cj^T^pLb) , 

{s-i^T''pLb)Y,{<il^Ty) , 
{sl^.,l^.,l^.,PLb)Y,i(il'''l^'l'''(l) > 



(2.23) 



e 

167r2 

9 
167r2 



sa>"'{mspL + mbpR)b 



F, 



where Pr^l = (1 =t 75)/2. The SUSY contributions to the b ^ sj transitions are encoded 
in the Wilson coefficients as additional contributions with respect to SM values. 

The values for the rare branching ratio Br[Bs — ?• in the SM to NLO are given 

by (see for example [28] and references therein) 



Br{Bs 
BriB, 



At > jSM 

e+e")sM 



3.1 X 10" 



2.4 X 10" 



w. 



ts\ 



X 



0.04^ 

BriB, 



0.21 GeV 



The current bound on the leptonic branching ratio is 



215 . 



^"^/J- jSM 



Br{Bs l-i^lJT) < 6.6 x 10"^ {2a) . 



(2.24) 
(2.25) 

(2.26) 



Note that the analogous i?^ channels will be suppressed relative to the Bg case by 
(/ly/IJ l^trf/^t^P ~ As we have already mentioned, in the NMSSM the light- 

est CP-odd boson ai can be as light as a few GeV and still satisfy constraints from 
Br{Bs — > fi'^fi"), for low to intermediate values of tan/3. 
We also check the contributions to the mass splittings 



AMx = 2IM12I 



1 + 



12 



M12 



(2.27) 
(2.28) 
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from B-mixing (AM^ and AM^). 

The isospin asymmetry Aq{B — t- K*j) is defined as 



7*0 



The ratio of leptonic kaon decays to pion decays in SUSY is 



(2.29) 



Br(K fiUf^) _ TK_ 
Br(7r iiu^) 



^1. 



rrid 



tan^ /3 
1 + e'g tan /3 



1 2 



(l + <5e 



(2.30) 



where is a loop factor and 6em is an electromagnetic correction term. R123 is defined 
as [35] 



123 



Vus{Ki2) Kd(0+ ^ 0^ 



K.(^Z3) Kd(vr«2) 



1 



m 



K+ 



1 _ "M^ tan^/? 



y 1 + Eg tan /3 



(2.31) 



where K12 and ET/s denote helicity-suppressed {K — ?• /z^) and helicity-allowed {K — t- vrZi/) 
decays. 0^ — ?• 0^ denotes nuclear beta-decay. Note that in the SM R123 = 1- 

There has been some discussion regarding Br{B — t- tv) and Br{B — t- Dtv) in the 
context of a relatively light charged Higgs and a possible tension between the current 
measurements of the two and with [g — 2)^ (For a recent paper discussing this aspect see 
[38]). The charged Higgs contribution in the MSSM may be written as [39, 40, 41] 



Br{B 



TV 



mlflms ( 1 



mz 



m 



B 



1 + 



m 



B 



NP 



where 



mbm£ tan^ /3 



m 



H+ 



1 + eo tan /3 



(2.32) 



(2.33) 



Rtut is the ratio between the measured Br{Bu — ?• TVr) and the SM prediction. The present 
Br{B — )• Tiy) measurements seem to suggest no/very small MSSM contribution. The SUSY 
contribution is generally substantial for large tan /3 or for light charged Higgs. 
The differential decay rate of i? — )• Divi can be expressed as [43] 



dViB ^ Dive) _ Gj,\Va\'m% 



dw 



1927r3 



m 



pv{w) 



1 - t{w) 



(2.34) 



B 



rub 



tan2/3 



(mb — mc)mjj+ 1 + eo tan /3 



ps{w) 



where pD, Pb are the meson four-momenta, w is defined as 



w 



1 + (mD/mB)'^ - (pB - PD)'^/m'^ 



B 



2m£)/mB 



(2.35) 
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Figure 2: In the MSSM, the charged-Higgs (H^) contributes to both B ^ tv and B Dtv 
in a similar way. The relevant Cabibbo-Kobayashi-Maskawa (CKM) factor in the former case is 
Vub which has a larger uncertainty than Vcb^ which is the relevant one for B — > Dtv. Apart from 
this B ^ TV involves two powers of the QCD form factor which must be obtained from lattice 
calculations and presently has a large uncertainty. The latter process involves two QCD form factors 
which are better constrained [42]. 

t{w) = m\ + m|, — 2w ttid uib and py{w), ps{w) are the vector and scalar Dalitz density 
distributions, respectively. 

The Br{B — )■ Dtv) may be parametrized using the same Wilson coefficient as [43] 



We would like to emphasize that the persistent indications from the muon anomalous 
magnetic moment discrepancy suggests some new additional contribution (maybe just from 
low-energy QCD, hadronic light-by-light or BSM, to speculate on a few sources). To 
achieve this in the purview of the MSSM/NMSSM one generally needs a sufficient tan/3 
(not necessarily very large though; we use tan /3 = 10, for instance, and still get sufficiently 
large contributions to {g — 2)^ in the Su interval). Hence we expect some tension between 
{g — 2)fj_ and some of the flavor observables in general. Thus imposing both of these will 
flush out the most interesting regions of the parameter space. 

2.3 Collider bounds 

The direct collider bounds come from various channels. For instance, a promptly decaying 



Br{B Dtv) 



(0.28 ± 0.02) [l + 1. 38(3) Re{CJ^p) + 0.88{2)\C]:^p\'^] . (2.36) 
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Lower bounds on the Higgs and sparticle masses (GeV) 




111 




73 


H+ 


79.3 




94 




93.4 


n 


81.9 


x\ 


46 


D 


94 


xl 


62.4 


ii 


95.7 


xl 


99.9 


h 


89 


xl 


116 


q 


379 


xt 


94 


9 


308 



Table 2: Mass bounds from collider studies [44]. 



NLSP places a limit on A = F/Mm in RGM of about 100 TeV. For non-prompt NLSPs 
the lower bound comes from charged sparticle masses and indirect constraints such as the 
inclusive tri-lepton signal. In both cases, but especially for promptly decaying NLSPs, a 
heavy spectrum results. 

The trilepton process can put indirect mass bounds. The process is 

PP ^ ixl ^ x?/+r ^ x'l^i^i) , (2.37) 

and proceeds mainly via an s-channel virtual W exchange (and a small t-channel squark 
exchange) and hence its cross section peaks at « Mw/2. Otherwise the cross section 
falls with increasing mass with a weak tan /3 dependence. This puts indirect mass limits 
on the masses. 

The current bounds on the superpartner masses [44] are shown in Table 2. For the 
sleptons, constraints from the decay width already put a lower bound of ~ 40 GeV 
independent of decay modes. The squark case includes cascade decays which gives a slightly 
lower bound than if the squarks were assumed to decay directly to photinos. Also, in 
deriving the squark mass bound with cascade decays, fixed values for fj, and tan/3 are 
taken. This is also true while deriving the g mass bound. It is important to note that these 
experimental bounds have been derived after imposing certain assumptions some of which 
are not strictly true in many studies, and definitely not in ours. Hence some of these must 
be viewed as indicative. Some of these assumptions include [44] : 

• Conservation of R-parity. 

• Gaugino mass unification at the GUT scale. 
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• xi is the LSP. 

• All scalar quarks are assumed to be degenerate in mass except for t and b. 

• rrig^ = ^QL 

• The mass of G, in the decay final states, is neglected relative to the other masses. 

• Some of the mixing angles are small. For instance, bounds from e"*" — e~ collisions 
depend on the mixing angle of the lightest squark mass eigenstate. It is assumed that 
only b and t have non-trivial mixing angles in this case. 

There are also model dependences, for instance, in NMSSM the mass of the lightest CP- 
odd Higgs can be quite small as opposed to a bound of 93.4 GeV. In Table 2 also note that 
we have taken the bound on the lightest CP-even Higgs, nih, to be 111 GeV rather than 
the usual LEP bound value of 114 GeV, to accommodate a possible ~ 3 GeV uncertainty 
in the theoretical Higgs mass calculation. 

3. Analysis 

3.1 Numerical analysis and strategy 

We have taken the mediation/messenger scale [Mm) to be high, in the vicinity of the GUT 
scale 

M„ ~ Agut . (3.1) 

The study may be extended for other messenger scales in a straightforward manner. Im- 
posing the messenger scale to be Aqut makes comparison with mSUGRA and gaugino 
mediation (with a compactification scale ~ Kqut) convenient and to lowest order removes 
any effects which are solely due to RG running from GUT scale to low scale (electroweak) . 
In this context the reader is also referred to [19, 20]. The characteristic scale associated 
with the hidden sector is Mg. With this definition the correlation functions {Ar^Bj.) be- 
come dimensionless. The GGM masses, Eq. (1.14), are imposed at the messenger scale 
as functions of the correlation functions {Ar^B^) and then RG evolved to the low scale. 
For completeness, we give the 1-loop RG equations for the MSSM and NMSSM in the 
Appendices. For the complete 2-loop RG equations the reader is directed, for example, to 
[45, 46]. 

One point we would also like to emphasize is that we do not take B^ as the fundamental 
quantity, in contrast to an interesting study in [8] for instance, and instead choose to 
keep tan /3 fixed at 10 to more readily compare with various benchmark scenarios in the 
literature. If Bji is taken as the fundamental quantity and to be vanishing at the GUT 
scale (by adopting a strict interpretation of gauge mediation) then tan/3 is no longer an 
input and comes out to be large (~ 15 — 65) as in [8]. A higher value of tan/3 can fit Aa^ 
for a heavier sparticle spectrum but for small to intermediate values for tan/? it is easier 
to accommodate flavor-physics constraints. We also adopt a pure gauge mediation setting 
where the assumption is that the generation of /i is independent of the SUSY-breaking 
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and the GGM scenario. As we previously mentioned the MSSM Higgs sector has 4 free 
parameters and the NMSSM Higgs sector has 7 free parameters. We use the minimization 
conditions to trade for tan /3. The // parameter is calculated at the low scale, using the 
minimization conditions, to be consistent with the input mass (M^o). Thus the actual 
inputs we choose in the MSSM Higgs sector are 

Mzo, tan 13,171^^,171^^ . (3.2) 

Here tan f5 is chosen at the scale of EWSB. Notice that Hu and if^ carry the same quantum 
numbers as l^. Thus we will assume that the predominant contribution to rn'jj^ and itt-h^ 
at the high scale are same as the II soft mass contributions. As previously mentioned 
benchmark value of tan /3 = 10 is adopted throughout our study, since our main focus is 
on the effects of the hidden-sector correlation function ratios. 

In the scale-invariant NMSSM case we take a slightly different set of inputs. The three 
Higgs potential minimization conditions in the NMSSM may be written as [21] 

Vd (^rnjj^ + filff + X'^vl + ^{vj- vl)^ - Vu f^eff. + '^^^eff) = , (3.3) 

Ml 

A 

Here 



Vu [mH^ + ^g^, + A + — (?;„ - v^)j - Vd fJ^ejj. [Ax + -fie 
m% + kA^^ + + >^'^ivl + vj) - 2XvuVd^ - XvuVdAx = . 



fieff. = A(iV) (3.4) 

The above minimization conditions may be used to eliminate three of the soft scalar mass 
terms. Using this freedom we choose the inputs 

M^o,tan/3,A,K,^A,^K,Me#. , (3.5) 

at the low scale. For our study we take tan/3 = 10 as before, A = 0.1, k = 0.1, A\ = 150, 
= and figff. is taken to be 210 GeV. Once again other suitable choice of parameters 
may be made, but our main focus in this study is on the dependence of the low-energy 
observables and spectra on the hidden-sector correlation function ratios. A very compre- 
hensive and detailed study of the NMSSM parameter space with various choices of the 
parameters may be found in [21, 47] and references therein. 

Since the MSSM contribution to the muon anomalous magnetic moment depends on 
sgn{fi), currently > seems to be slightly favored and we focus on a positive value 
throughout. 

We have used the following codes in their original and modified forms : SuSpect [48], 
Superlso [36, 49] and micrOMEGAS [52, 51, 50]. We have cross-checked the theoretical 
predictions of the low-energy observables between the programs and also independently 
estimated some of their values directly from the low-energy spectra. For the NMSSM we 
use our code to implement GGM soft masses at the GUT scale and evolve it down to a low 
scale (Electro- weak) . The NMSSM RGEs are implemented only to 1-loop (see Appendix 
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B.) and are found to introduce an error in the spectra of atmost ~ 10 — 15%, mainly in 
the third generation, in comparison with mSUGRA 2-loop NMSSM implementations (for 
example in NMSSMTools [53, 54]). At the low scale we have also used NMSSMTools [53, 54] 
to calculate and cross-check some of the low-energy observables. 

3.2 Features of the parameter space 

We parametrize the various correlation function ratios with the notation 

(^1 : B2 : Bs \ A, : A2 : A3) , (3.6) 

where the and A^ are defined as in Eq. 1.14. There are a couple of caveats. Note 
that there could be a hierarchy between the B,. and Ar correlation functions even when 
the case under consideration is (1 : 1 : 1|1 : 1 : 1). We elaborate on this below. We also 
remind ourselves that in order to make the Br and A^ dimensionless, we have factored out 
a common scale Ms = 2x 10^ GeV, that may be thought of as a characteristic scale for 
the hidden sector. But this is just our convention. 

To keep the investigation of low-energy observables manageable we will mainly focus 
on scenarios where the ratios are mirrored between the Br and Ar ( i.e. the ratio associated 
with the SU{2)l and 5?7(3)c correlation functions with respect to the U{1)y correlation 
functions are the same for both Br and Ar). Other scenarios may be considered in a 
straightforward manner. The viable GGM parameter space for three correlator ratios 
(1 : 1 : 1|1 : 1 : 1), (1 : 3 : 3|1 : 3 : 3) and (1 : 1/3 : 1/3|1 : 1/3 : 1/3) are shown in Figs. 3 
and 4. As previously mentioned, the RGM limit can be obtained from the (1 : 1 : Ijl : 1 : 1) 
GGM scenario by considering points where 




At these points, in the (1 : 1 : 1|1 : 1 : 1) case, there is no hierarchy between the Br and A^. 
Large parts of the GGM parameter space are excluded purely from requirements of having 
no tachyonic masses, potential to be bounded from below, absence of charge breaking and 
perturbative RG evolution. Nevertheless there are viable regions different from RGM that 
have interesting observables and spectra. 

There are also non-viable regions in Figs. 3 and 4 where the value of the MSSM 
contribution to {g — 2)^ is above 3ct relative to the current discrepancy between theory 
and observation. It is also found that large regions are excluded due to flavor-physics 
observables with a large overlap with exclusions due to direct collider bounds. Direct 
collider searches, flavor-physics and muon anomalous magnetic moment give interesting 
bounds on the gaugino and sfermion masses in the GGM scenario. 

Larger gaugino masses can cause the sfermion masses to run to larger values due to 
RGE and hence lead to more phenomenologically viable regions. Looking at the plots, it is 
interesting that log^Q Ar does not seem to be bounded from below. This may be understood 
by noting that, for fixed gaugino soft mass terms, the ^4^ — >• limit implies that at Agut 
we are basically setting all m? — )■ (since we have set the Fayet-Iliopoulos term C = 0, by 
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Allowed GGM parameter space (MSSM) 



B„. 1:1:1, A . 1 :1 :1, Tan p= 10, M = A, 



J ■ 



Parameter space 
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[A a f'^" 2 3a [A a f^' 


■ 


Tachyonic/ NaN mass 




RGE NP/Landau pole 




Direct searches 


■ 


flavor-exclusion 






Log 



Allowed GGM parameter space (MSSM) 




Log 

Figure 3: Constraints on the (1 : 1 : 1|1 : 1 : 1) (top) and (1 : 3 : 3|1 : 3 : 3) (bottom) GGM 
parameter spaces (MSSM). Points in the vicinity of the Hnear band (black) for (1 : 1 : 1|1 : 1 : 1), 



where log^o ^> 



(1/2) log;^o^a' correspond to RGM scenarios. For the (1 : 3 : 3|1 : 3 : 3) 



case, regions near the hnear band have no RGM interpretation and the band is merely given for 
comparison. In the figure the acronym RGE NP stands for non-perturbative renormalization group 
evolution. The exclusions due to {g — 2)^ are shown in red, those due to flavor observables in 
magenta/pink and those due to direct collider bounds in grey. The white portions represent the 
allowed regions. 



say messenger parity arguments). The initial RG running of the sfermion masses from the 
GUT scale are therefore completely determined by the gaugino soft terms (see Appendix 
A) and hence the sfermion masses at the low scale approach a constant value as we decrease 
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Figure 4: Constraints on the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) GGM parameter space (MSSM). Note 
again that for this ratio, regions near the Hnear band (black) do not correspond to RGM. The white 
portion again represents the viable region. 



In RGM there is very little freedom to change the sfermion to gaugino mass ratio 
except by 0{1) factors. The ratio may be tuned by varying the number of messenger 
fields (non-minimal RGM), but larger ratios may be disfavored since they may spoil gauge 
coupling unification. In GGM we can explore larger ratios between them generically and 
even among the three gauge groups, not restricted just by ratios of g'^, due to the possibility 
of hierarchy among the three Br ■ 

Let us take up the question of precision gauge coupling unification in the MSSM. There 
is a mismatch in of about 3% at the GUT scale derived from the unification of ai and 
a2- The RG evolution of the coupling constants, taking into account threshold effects, obey 

dgi hi 



dt 167r2®^^^'' (^i'&2,M = (f,l,-3) 



(3.8) 



with t = ln(g//i). The 0j take into account the thresholds due to the superpartner masses 
and are given by 



1 

33 



i=l 



10 

y 



3 



(3.9a) 
(3.9b) 



=1 



3 3 18 ^ V 

i=l 



where the Heaviside functions above are defined as 
all the thresholds, the 0j — >• 1 (see Appendix A ). 



(3.9c) 

- mp . At a scale much above 
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fx'-(g-2) constraints on GGM in the MSSM 




^x^-(g-2) constraints on GGim in tlie imSSIVI 




Figure 5: Points that give a positive contribution to {g — 2)^, for allowed regions in the (1:1: 
1|1 : 1 : 1) (top) and (1 : 3 : 3|1 : 3 : 3) (bottom) cases (MSSM). Points in the vicinity of the linear 
band (black) for (1:1:1|1:1:1) correspond to RGM as before. The color scale to the right shows 
the MSSM contribution to {g — 2)^ in units of 10^^^. Note that the color scales are different in 
each case. The point marked with the star symbol is an example point that we will explore later, 
as a benchmark, in the context of low-energy observables and mass spectra. 

A straightforward examination of the above expressions tells us that a heavy wino (w) 
and/or a light gluino (g) can lead to precision gauge coupling unification by changing the 
running of a2 and/or 03. The other colored particles can also slow down the running of 03 
but they come with a smaller coefficient and moreover will also lead to changes in the ai 
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and 02 running which makes it more tricky to implement. It is more efficient to work with 
a hght g as opposed to a heavier w to achieve precision gauge couphng unification since 
this way one can utiHze the larger value of 03 at lower energies. 



^x'-(g-2) constraints on GGM in the MSSM 




Figure 6: MSSM contribution to [g ~ 2)^ for the allowed region in the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) 
case. The color scale to the right show the value of (5— 2)^ again in units of 10^^^. The regions ruled 
out by collider bounds are shown in grey, those ruled out by flavor physics by pink/magenta and 
those with a (g — 2)^; > 3ct by red. As in the (1:3: 3|1 : 3 : 3) case the linear black band is merely 
given for comparison with the (1;1:1|1:1:1) case and does not have an RGM interpretation. 
The star symbol again denotes the benchmark point. It is ruled out in the present case from collider 
and low-scale observable bounds. In fact, note that for the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) case regions 
below logio-^i ^ are ruled out in contrast to the (1 : 1 : 1|1 : 1 : 1) and (1 : 3 : 3|1 : 3 : 3) 

cases. 

To lowest order the achievement of precision gauge coupling unification therefore cru- 
cially depends on the mass ratios of the gauginos (hence on i?i : B2 ■ B3), specifically 
that of the w and g (note that the bino {B) does not contribute to the running). Thus 
in the GGM scenario there is more room for variations in these ratios to possibly improve 
precision gauge coupling unification at the GUT scale. 

Figures 5 and 6 show the contribution to the muon anomalous magnetic moment 
for the same correlation function ratios (1 : 1 : 1|1 : 1 : 1), (1 : 3 : 3|1 : 3 : 3) and 
(1 : 1/3 : 1/3|1 : 1/3 : 1/3). Note that {g — 2)^ places strong constraints on the space if it is 
assumed that the MSSM contributions are viable explanations for the current discrepancy 
(i.e. the contribution is within at least 3a of the current discrepancy between SM theory 
and experiment). 

In the (1 : 1 : 1|1 : 1 : 1) case the region in the interval logj^g -^1 ^ [—6.5,-6.3] 
and log]^o ^1 ^ —13.5] have {g — 2)^ contributions that are in the interesting range 

for explaining the current discrepancy. Regions with log^g Bi < —6.5 have even larger 
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u -(g-2) vs. b ^ sy 




^1 -(g-2) vs. b ^ sy 
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Figure 7: Aa^ vs. BR(6 ^ sj) for the (1 : 1 : 1|1 : 1 : 1) (top), (1 : 1/3 : 1/3|1 : 1/3 : 1/3) (right) 
and (1 : 3 : 3|1 : 3 : 3) (left) cases (MSSM). Note again that the direct colhder bounds have not 
been imposed. We note that there is a slight anti-correlation between the {g — 2)^; contribution and 
the BR(& — >• s'-f), which is most obvious in the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) case. 

{g — 2)^ values but are ruled out by collider and flavor-physics bounds. In regions outside 
logxo -Bi ^ [—6.5,-6.3] and log^o e [-18,-13.5] the {g — 2)^ contribution drops rapidly 
and have low values that are not sufficient to explain the discrepancy. 
Changing the correlation function ratios from 

(1:1: 1|1 : 1 : 1)mssm — ?" (1 : 3 : 3|1 : 3 : 3)mssm j 

drastically reduces the {g—2)f^ values in the log^g Bi e [—6.5, —6.3] and log^g ^ [~18, —13.5] 
intervals. The values are now lower by a factor of 3-6 and unlikely to explain the discrep- 
ancy by themselves. 

For the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) case regions below log^g i?i < —6.1 become non- 
viable due to {g — 2)^ > 3a, flavor-physics exclusions and collider bounds. The interesting 
region is now shifted slightly to log^g -^i ^ [~6.1, —5.9] and log;^o ^ [~18, —13.0]. 

Thus in the MSSM a combination of low-energy observables and consistency conditions 
place severe restrictions on the GGM parameter space. We will later look at the NMSSM 
case and come to a similar conclusion. In conclusion, for the (1:1:1|1:1:1) and 
(1 : 1/3 : 1/3|1 : 1/3 : 1/3) cases, there is a very narrow band where the {g—2)^, contribution 
is within the 2(j interval of the current discrepancy. 
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NLSPs in MSSM/GGM 
Gluino NLSP 

■ Sneutrino NLSP 

■ Slepton NLSP 
Squark NLSP 
Neutralino NLSP 
Chargino NLSP 




Figure 8: NLSP species in the (1 : 1 : 1|1 : 1 : 1) (top) and (1 : 3 : 3|1 : 3 : 3) (bottom) GGM cases 
(MSSM). It is interesting tliat the NLSP topography looks very similar in the (1 : 1 : 1|1 : 1 : 1) 
and (1 :3:3|1 :3:3) cases in spite of differences in the over-all mass spectra and other low-energy 
observables. The benchmark point we take as an example is again shown by the star symbol. 



The Aa^ vs. Br{B — t- ^^7) plot is shown in Fig. 7, for the three correlation function 
ratios. It is noticed that in the (1 : 1 : 1|1 : 1 : 1) and (1 : 3 : 3|1 : 3 : 3) cases 
the Br{B — t- X^^) is confined to a thin sliver between the values [3.2,3.4] x 10~^. For 
intermediate values of {g — 2)^ the allowed Br{B — t- Xs'j) interval is the broadest. In 
the (1 : 1/3 : 1/3 11 : 1/3 : 1/3) case the Br{B — t- Xsj) interval enlarges appreciably to 
[2.0, 3.4] X 10~^ for allowed values of {g — 2)^. We note a slight anti-correlation between Aa^ 
and Br{B — ?■ ^^7) in this case. There is also an aggregation of points for Br{B — )• X<j7) 
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in the interval [3.2,3.4] x 10"^. 



NLSP regions (MSSM) 
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Figure 9: NLSP species in the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) case (MSSM). The NLSP topography 
is quite distinct from the (1 : 1 : 1|1 : 1 : 1) and (1:3: 3|1 : 3 : 3) cases. The collider exclusions 
are shown in grey and the flavor exclusions in magenta/pink. The flavor excluded regions have an 
overlap with regions excluded by {g — 2)^ > 3a and are not visible in the figure. 

Let us now consider the NLSP species in each of the cases. Figs. 8 and 9 survey the 
topography of the NLSP species. In the (1 : 1 : 1|1 : 1 : 1)mssm and (1 : 3 : 3|1 : 3 : 3)mssm 
cases the most favored NLSP candidate in the region with interesting values of {g — 2)^ 
(i.e. within 2a of the current discrepancy) is the stau (f) with a small region (around 
logj^o -^1 ~ [—6.5,-6.3], logxo ^1 ~ [—15,-13]) that also accommodates neutralinos (x^). 
There are also regions with small contributions to (g — 2)^ that have both the slepton and 
neutralino as the NLSP. Observe that in the RGM case the NLSP would have been the 
solely and would give a very low value for {g — 2)^ in that region. In this context also note 
that there are viable regions with Mx <C mj with a x^ NLSP which nevertheless give a very 
low contribution to {g — 2)^. Do note that for the above two cases there are boundaries 
where the xp and / are almost degenerate in mass. In these regions x^-l coannihilations 
might be important [5] and could lead to interesting phenomenology. Also note that there 
are no viable regions with a chargino {x^) or sneutrino (P) NLSP in these two cases. 

In the (1 : 1/3 : 1/3)1 : 1/3 : 1/3)mssm case the NLSP topography is completely 
transformed. In contrast to the (1 : 1 : 1)1 : 1 : 1)mssm and (I : 3 : 3)1 : 3 : 3)mssm 
cases there are viable regions with a x^ or u NLSP. In the regions with a significant 
{g — 2)fj_ contribution, the favored NLSP is now solely the x^. In the approximate interval 
logio^i e [-5.2,-4.3] and log^o -4i e [-18,-11.0] the D is the NLSP. For large values of 
logj^o ^1 ~ — ^5 between the intervals log;^Q Bi ~ [—5.5, —4.1] the x^ along with the x^ are 
degenerate in mass. Thus in contrast to (1 : 1 : 1|1 : 1 : 1)mssm and (1 : 3 : 3|1 : 3 : 3)mssm 
now there are regions and boundary regions with x^ — x^ and — i' degenerate in mass. 
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Figure 10: Aa^ vs. TOnlsp for the (1 : 1 : 1|1 : 1 : 1) (top), (1 : 3 : 3|1 : 3 : 3) (middle) 
and(l : 1/3 : 1|1 : 1/3 : 1/3) (bottom) cases (MSSM). As previously mentioned we only focus on 
positive {g — 2)^ contributions. Note that the direct collider bounds have not been imposed in the 
plot. The constraints from the collider bounds, in each of these cases, can nevertheless be easily 
read off by comparison to the (g — 2)^ plots presented earlier. The filled regions denote allowed 
values. 



In the vicinity of these regions again co-annihilations among the almost degenerate NLSP 
candidates become important and will lead to distinct phenomenology. 

A plot of Aa^ vs. mNLSP is shown in Fig. 10 for the three correlation function ratios. 



The filled regions denote the allowed values. It is clear that for each of the cases the allowed 
values of the NLSP mass are correlated, albeit weakly, with the {g — 2)^ value. We observe 
that for large values of the {g — 2)^ contribution the allowed range of the NLSP masses is 
smaller than for lower values of {g — 2)^. In the (1 : 3 : 3|1 : 3 : 3) case the NLSP mass is 
within a smah interval of [100,275] GeV. For the (1 : 1 : 1|1 : 1 : 1) and (1 : 1/3 : 1/3|1 : 
1/3 : 1/3) cases the NLSP mass is within a broader interval stretching all the way from 
100 GeV to ~ 600 GeV. For a given value of {g - 2)^ the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) case 
accommodates a wider range of allowed NLSP masses than the (1 : 1 : 1|1 : 1 : 1) case. 

The mass of the lightest CP-even Higgs is plotted in Figs. 11 and 12. For the (1:1: 
1|1 : 1 : 1), (1 : 3 : 3|1 : 3 : 3) and (1 : 1/3 : 1/3|1 : 1/3 : 1/3) cases in the region with 
interesting values of {g — 2)^, the mass of the lightest CP-even Higgs (including the leading 
2-loop corrections) is found to be typically in the interval [114, 125] GeV. The actual value 
at a given point in parameter space is seen to depend sensitively on the correlator ratios. It 
is found in general that increasing the ratio of the other correlation functions with respect 
to Bi and Ai increases the Higgs boson mass and conversely decreasing the ratio decreases 
the mass. 

RGM models with very large Nm have a superpartner spectrum similar to gaugino 
mediation (where the scale is typically taken close to the GUT scale). GGM can therefore 
be thought of as an interpolation between RGM and gaugino mediation [5]. This is also 
another motivation for taking our messenger scale to be near Agut, so as to make it more 
convenient to compare GGM scenarios with mSUGRA or gaugino mediation scenarios. 

Owing to the freedom of current- current correlator ratios there may be scenarios where 
GGM gives a squeezed superpartner spectrum alleviating the little hierarchy problem. If 
is lower with respect to B2, for example, this will compress the low-energy spectrum (mak- 
ing for instance the g-x mass ratios smaller) while if it is higher it will further increase the 
hierarchy. For example, at the point (—6.5, —17.5) (MSSM case) in the (log^g -^i) loSio ^1) 
space (with tan/3 = 10 and Mm = Aqut) going from 

(1:3: 3|1 : 3 : 3)mssm — >■ (1 : 3 : 1|1 : 3 : 3)mssm j 

changes the g-x mass ratio by 




2^ 1, 
16 ^ 6 . 



Lowering B2 with respect to B^ can similarly make the mass ratio selectively larger. 

The RGM intuition of large squark {q) to slepton [l) mass ratios need not again nec- 
essarily hold true in GGM. This is because generically the I have larger hypercharges and 
for large enough Ai the hypercharge term may selectively give a larger contribution to I 
masses as compared to the squarks. Going from 

(1 : 1 : 1|1 : 1 : 1)mssm (1 : 1 : 1|/C : K. : /C)mssm , 
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Figure 11: Mass of the lightest CP-even Higgs in the MSSM for the (1 : 1 : 1|1 : 1 : 1) (top) and 
(1 : 3 : 3|1 : 3 : 3) (bottom) GGM cases. The color scales are in units of GeV and different in each 
case. 



at the high scale with /C ^ 1 is seen to decrease the q to / low-scale mass ratio appreciably 
in large parts of the viable GGM space. Note that this is equivalent to raising the over-all 
scale of the Ar correlation functions relative to the scale of the correlation functions. 
For example starting at (—6.0, —12.0) (which for (1 : 1 : 1|1 : 1 : 1) is an RGM like point) 
in the (log^o -^i' loSio ^i) space, with tan/3 = 10 and Mm = ^gut as before, going from 
(1 : 1 : 1|1 : 1 : 1)mssm to (1 : 1 : 1|/C : /C : /C)mssm with /C = 20 changes the q to / ratios at 
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Figure 12: Mass of the lightest CP-even Higgs in the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) case. Note that the 
color scales are again in units of GeV and different from the (1:1:1|1:1:1) and (1:3:3|1:3:3) 
cases. Once again the collider bounds are in grey, the flavor bounds in magenta/pink and {g — 2)^ 
exclusions in deep red (unhatched). 



the low-energy scale from 



ma , 

" - (3.10) 



Existing models of gaugino-mediation have large hierarchy between sfermion and gaug- 
ino masses at the compactification scale (typically ~ Aqut) which gets washed out at EW 
scales due to RG running. For instance a typical spectrum in gaugino mediation at the 
compactification scale (~ ^gut) niay be (see for example [14, 56]) 

Ma = mi/2 Vael,2,3 , 

i2 



m ? 



2 ^ '"1/2 



/ ~ 167r2 ' 



A 



167r2 ' 



which exhibits a huge gaugino to sfermion mass ratio of C'(47r). This ratio gets reduced to 
0(1) when we run the masses from Kqut to Kew [56] . If the compactification scale is 
taken to be low, then such a "low-scale" gaugino-mediation usually prefers a I NLSP over 
a B NLSP. This is due to the short RG running scale. 
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Figure 13: A section of the viable parameter space in the NMSSM for the (1 : 1 : 1|1 : 1 : 1) case 
showing the Aa^ contribution (top) and the NLSP candidates (bottom). Regions with tachyonic 
masses or where the sohitions to the 1-loop RG equations failed to converge are shown in dark 
grey. Light grey represents regions excluded by collider bounds, brown denotes regions where 
(g — 2)^ is above 3a and magenta represents regions ruled out by flavor-physics constraints. The 
Aa^ contribution for this choice of parameters is relatively small. Note that the lightest CP-odd 
scalar in this case generally has a low mass (see for example the benchmark point in Table 4) and 
its 1-loop contribution is negative. The regions near the linear black band correspond to RGM and 
it is clear from the plots that, with the current choice of model parameters, most regions with an 
RGM-like spectra are already excluded. The NLSP is the neutralino (represented in green) in the 
allowed regions shown. 

Generically, GGM may favor non-split {M\ ~ m^) or mildly split SUSY (in split SUSY 
we have Mx ^ characteristic of RGM and direct/hybrid gauge mediation models 

respectively. For instance as pointed out already in [5] and is clear from our figures there 
are points in the GGM space that have a spectrum with gauginos lighter than sfermions. 
These regions in contrast to RGM could give rise to very different phenomenologies. This 



-28- 




-8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 



Log 

Figure 14: The lightest CP-even Higgs mass in the (1 : 1 : 1|1 : 1 : 1) case. The mass is relatively 
small and in the [111, 125] GeV window. As in the previous figure, regions with tachyonic masses 
or where the solutions to the 1-loop RG equations failed to converge are shown in dark-grey. Light- 
grey represents regions excluded by collider bounds, brown denotes regions where (g — 2)^ is above 
3(7 and magenta represents regions ruled out by flavor-physics constraints. 

is characteristic of 'direct gauge mediation' models (where messengers participate in SUSY 
breaking) which usually have gaugino masses suppressed with respect to scalar masses. 

In Figs. 13 and 15 we show for the NMSSM the distribution of (g - 2)^, NLSP and 
lightest CP-even Higgs mass in a section of the viable GGM parameter space. For our 
choice of parameters it is observed that the {g — 2)^ is not very large and in the viable 
regions the is the NLSP. 

In Table 3 we show the dependence on the correlation function ratios for a benchmark 
point, in the MSSM, chosen so that the lightest CP-even Higgs mass is near the LEP 
bound. It is noticed that for this benchmark point in going from (1:1: 1|1 : 1 : 1)mssm 
(1 : 3 : 3|l : 3 : 3) MSSM the flavor-physics contributions are unchanged but the value of 
the contribution drops roughly by a factor of 4. This may be understood partially by 
noting that the x^, x^ masses have increased. Thus as alluded to, while discussing the 
similar NLSP topography in the two cases, the low-energy observables are quite different 
at each point in the parameter space. Table 4 shows similarly a benchmark point for the 
NMSSM case. The choice of the benchmark points in both the MSSM and NMSSM were 
guided mainly by requirements of a viable mass spectra and low-energy observables, but 
we have specifically picked a point which gives a relatively light CP-even Higgs (close to 
the current LEP bound). 



Regions where mj <^ M\ may imply that the Poppitz-Trivedi type mass terms (these 
are R-symmetry preserving terms accounted for by a non- vanishing messenger supertrace) 
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Description/Units 


Quantity 


(1 : 1 : 1|1 : 


1:1) 


(1 : 3 : 3|1 : 3 


:3) 


Inputs (High-scale): 


logio 


-6.4 




-6.4 






logio ^1 


-17.5 




-17.5 








^GUT 








1 n T^T 1 -j-c; / 1 ,r\\lCT Gr* Q \pi\ 


fan R 


10 


10 


Low energy 


Attn 


139 X 10" 


-11 


34.9 X 10" 


11 


contributions: 


BR{B Xsj) 


3.02 X 10 


-4 


3.10 X 10- 


-4 




BR(B, a^uT) 


3.4 X 10- 


-9 


3.16 X 10- 


-9 




\ (Jj 1 / 


9.3 X 10- 


-5 


9.41 X 10- 


-5 




BR{B D^TVr) 


6.9 X 10- 


-3 


6.91 X 10- 


-3 


giuino (^vjC V j . 




945.3 


2608 


charginos: 


m ± 

Xi 


308.2 




980.7 




(GeV) 


X2 


548.7 




1399.4 




neutralinos: 


mo 

Xi 


164.2 




167.1 




(GeV) 


m^o 


308.4 




980.7 






in 


533.2 




1391 






TM n 

A4 










Higgs: 


m,jO 


114.3 




121 




(GeV) 


m^o 


588.6 




1578.8 






m^o 


588.3 




1578.7 






m^^i 


594 




1581 




squarks: 




654.3 




1841 




(GeV) 




848.9 




2204.7 








862.5 




2355 








833 




2243 






mr 


866 




2356.2 






mr 

<1r 


830.6 




2244 




sleptons: 




147.4 




128.5 




(GeV) 


mf2 


276.8 




763.7 








274.3 




764.6 






mfiR 


155.7 




149.8 








263.1 




760.7 








262.3 




759.1 





Table 3: A benchmark point in tlie GGM parameter space for two correlation function ratios 
(MSSM). The point was chosen so as to have a light Higgs very near the LEP limit in the (1:1: 
1|1 : 1 : 1) case. This point is completely ruled out in the (1 : 1/3 : 1/3|1 : 1/3 : 1/3) case from 
collider, flavor and {g — 2)^ bounds. In both the cases above, f is the NLSP. Among other things, 
also note that changing the ratio from (I : I : I|I : I : I) — >■ (I : 3 : 3|I : 3 : 3) causes the {g — 2)^ 
contribution to decrease by a factor of 4, with a corresponding increase in the masses. 
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Description/ Units 


(ojuantity 


(1:1:11:1:1) 


(1:2:21:2:2) 


Inputs (High-scale): 


logio ^1 


-6.8 


-6.8 




1 A 

logio ^1 


-12.1 


-12.1 




Mm 


A 

J^GUT 




Inputs (Low-scale) 


A 


0.1 


0.1 






U.l 


0.1 






150 


150 














/^e# (GeV) 


zlU 


01 n 




-4- n /3 

tan p 


1 n 
iU 




IjOW energy 


An 


Al V 10-11 


97 V 1 n— 11 


contributions : 


BR(B X.'j) 


3.34 X 10""^ 


3.45 X 10""^ 






3.51 X 10"^ 


3.53 X 10"^ 




BR(5„ ^ rz/^) 


13.1 X 10"^ 


13.1 X 10"^ 


gluino (GeV): 


M3 


542.8 


1134 


charginos: 


m ± 


113.8 


191.1 


(GeV) 


m ± 
X2 


257 


353.3 


neutralinos: 


Xi 


62 


38.7 


(GeV) 


X2 


116 


192.5 






223.7 


223.4 






256 


353 




m^o 

A 5 


428.0 


432 


Hiffgs: 


777,1,0 


114.1 


117 


fGeV) 


777 1,0 


420.2 


420.3 






853.7 


846.4 






10 7 
iU. 1 


Q 




000.0 






777^± 


857.3 


850 


squarks: 




iUo4 


zUoU.o 


(GeVj 




iouy. ( 


Olio 
ZiiZ 


q.t. 




ioyi.o 


ioyo 






izoi. / 








ioyo.o 


1 QO/I 

ioy4 




777T 
<1r. 


1170. z 


1562 


sleptons: 


777fi 


910. 


DOO.O 


(GeV) 


777f2 


915.5 


935.7 






915.2 


1106.2 






916 


1128 






912.7 


1125.5 






911 


932.5 



Table 4: Benchmark points for two correlation function ratios (NMSSM). The point was again 

chosen so as to have the mass of the lightest CP-even Higgs boson near the LEP limit. For larger 
correlation function ratios, for instance (1 : 3 : 3|1 : 3 : 3) or (1 : 1/3 : 1/3|1 : 1/3 : 1/3), this point 
is again found to be non-viable. Note that we have only implemented 1-loop RGEs (Appendix B ) 
for this estimation as mentioned before. 
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Log 




Log B, 

Figure 15: A section of the viable parameter space in the NMSSM for the (1:3: 3|1 : 3 : 3) 
case, showing the Aa^ contribution (top) and the NLSP candidate (bottom). The color schemes 
are identical to the (1 : 1 : 1|1 : 1 : 1)nmssm case. Changing the correlation function ratios from 
(1 : 1 : 1|1 : 1 : 1) to (1 : 3 : 3|1 : 3 : 3) has further decreased the {g — 2)^ contribution but not 
appreciably. Also note that certain regions previously allowed in the (1 : 1 : 1|1 : 1 : 1) case are 
now non-viable and there are new regions that are viable. The NLSP species in this section of the 
viable parameter space is still the neutralino. 

in these regions are small [4]. In our specific context the Mx < may be understood 
as a consequence of a hierarchy between the Br and A,, correlation functions, for all other 
parameters remaining the same. Thus by raising the overall scale of the Ar relative to 
the Bj. the spectra may be inverted from an RGM like spectra. For example, consider the 
starting point (—6.5,-13.0) in the (log^g -^i) log^o ^1) space (again with tan/3 = 10 and 
Mm = ^gut) which would roughly correspond to RGM in the (1 : 1 : 1|1 : 1 : 1)mssm case. 
Then if we consider the re-scaling of the A^. correlation functions as 

(1 : 1 : 1|1 : 1 : 1)mssm — ^ (1 : 1 : 1| /C : /C : /C)mssm , 
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Log 

Figure 16: The lightest CP-even Higgs mass in the (1 : 3 : 3|1 : 3 : 3) case. 

for /C ^ 1 as before, this changes the gluino-sfermion mass ratios to give an inverted 
spectrum 




< 1 . 



There also exist even at the low scale other very approximate scaling relations, between 
the superpartner masses and the hidden sector correlation function ratios, in large parts 
of the viable parameter space that are not completely washed out by the RG running. For 
instance it is found that in regions with small values of Ar for 

(1 : 1 : 1| 1 : 1 : 1) ^ (/C : /C : /C| 1 : 1 : 1) , (3.11) 

at the high scale, we have the approximate scaling 

{mg,rn^Q,m^~j^,mq,mj,mi^)j^^;v ^ ~ (/C m^, /C m^-p, /C m^-^ , /C m-g, /C /C mc;)Ew • (3-12) 

This may be understood by noting that when the Ar are relatively small the initial RG run- 
ning of the sfermion masses is dominated by the gaugino mass terms (see Appendix A and 
B). The above scaling is just tantamount to changing the over-all scale of the Br correla- 
tion functions, and hence the gaugino masses at the high-scale. To give a concrete example 
consider in the (1 : 1 : 1| 1 : 1 : 1) case the point (—6.0,-16.0) in the (log^g log^Q Ai) 
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space. Then we have the following at the low scale : 

m-g (GeV) : 2215 'Szlll^ nge 4223 9951 

m^o (GeV) : 434 209 898 2339 , 

m^o (GeV) : 1208 'SzlII^ 665 2206 4903 , 

^nghtcBt ^Qgy^ . g-^g ^^=1/2^ ^ggg 4226 , 



'C=l/2 /C=2 ^^^^ /C=5 



wvicst /QgyN . -^207 ^ 665 2206 4903 , 

A 



lightest 

m- 


(GeV) : 


1570 


heaviest 
'"-(J 


(GeV) : 


2007 


lightest 

rrij 


(GeV) : 


367 - 


heaviest 
'"'I 


(GeV) : 


667 - 


lightest 

m- 


(GeV) : 


662 - 


heaviest 


(GeV) : 


663 - 



5zlII^ 815 ^ 2991 ^ 6983 , 
^zlZ^ 1065 ^ 3794 ^ 8835 



^zlZ?^ 184 ^ 729 ^ 1809 , 
342 ^ 1312 ^ 3207 





> 




=2 




— > 




=2 








=2 
— > 



^^-^ 331 1306 3198 





— > 


IC- 


=5 




— ;> 




=5 
— !> 



^^-^ 332 1310 3205 . 

To conclude, in this section we explored the effects of different hidden-sector correlation 
function ratios on low-energy observables and the SUSY mass spectra. Our main guiding 
points were observables from flavor physics and the muon anomalous magnetic moment. 
It was found that large regions of the GGM parameter space were disfavored (though not 
necessarily ruled out) from the viewpoint of these and there were interesting regions that 
gave sizeable contributions to {g — 2)^ on par with the current discrepancy. 



4. Summary 

It is indicative from the present study and some of the earlier ones [8] that low-energy 
observables can play a very complementary role to collider studies in deducing viable and 
interesting regions in the GGM parameter space. Specifically, we saw that there are strong 
relationships between the correlation function ratios in the hidden sector and the values of 
the low-energy observables and mass spectra at the low scale. Let us recapitulate some of 
the salient features noticed in the GGM case: 

• Muon anomalous moment and flavor physics place strong constraints. The detailed 
features of the parameter space are found to depend sensitively on the correlation 
function ratios. 

• The GGM scenario can more easily accommodate precision gauge coupling unifica- 
tion. This is due to the larger freedom of non-universal gaugino masses (due to 
different Bj.) that could lead to a light g or heavy w relative to RGM scenarios. 
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• The topography of NLSP is found to depend on the correlation function ratios. Even 
in cases where the NLSP topography is naively similar there are significant differences 
in the mass spectra and low-energy observable values. A case in point is the (1:1: 
1|1 : 1 : 1) and (1 : 3 : 3|1 : 3 : 3) cases that we considered. Also, it is found that there 
are interesting regions and boundaries in the parameter space with multiple NLSPs. 
In these regions co-annihilations may be important and could lead to interesting 
phenomenology. 

• We find that in most regions with interesting values of {g — 2)^ (i.e. within 2a of the 
current discrepancy) the mass of the lightest CP-even Higgs boson is relatively small. 

• There are very approximate, but interesting, scaling relations observed in large parts 
of the GGM parameter space (between the correlation function ratios at the high 
scale and sparticle masses at the low scale) that are not completely washed out by 
RG running. 

• There are viable regions where the gauginos arc lighter than sfcrmions. These regions 
may have interesting phenomenological implications. This may be understood as a 
hierarchy between the Br and Ar correlation functions. 

• There are also allowed regions where the RGM intuition of large q to I ratios is no 
longer true. These are regions where the correlation functions Ar are large. These 
regions favor larger I masses relative to q due to the fact that I hypercharges are 
typically larger than the corresponding q hypercharges. Thus there is a slight en- 
hancement to the I masses in these regions. 

• The ^ to X mass ratio may again be small in contrast to RGM expectations. This 
may be partially understood as a consequence of the hierarchy between B2 and 
correlation functions. 

• The Ar arc unbounded from below. This is the limit where the sfermion soft masses 
at the GUT scale tend to zero, since we have put the Fayet-Iliopoulos term C = 0, and 
the initial sfermion mass generation is essentially due to the gaugino masses during 
RGE. 
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Appendix A. 1-Loop MSSM renormalisation group equations 

For completeness we list the 1-loop RG equations for the MSSM and NMSSM. The RGEs 
are written in the third family dominant approximation with the definitions t = log{Q^/ Qq) 
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and U{1)y gauge coupling constant gf = ^{gi^^)"^- The terms in the MSSM superpotential 
and the soft terms are as defined in Eqs. (1.2) and (1.3). For the complete 2-loop RG 
equations, see for instance [45, 46]. 



A.l Gauge and Yukawa couplings 



167r 



167r 



IGvr- 



2 dgj 
dt 
2 ^91 
dt 

dt 



9l 



(1) 



dt 

A. 2 Gaugino masses 



16-^ = (6^' + y" +yr- l^l - ^^l - "^gl 

W7r'^ = y'My'r + 3y!-3gl-3gl) (2) 



IGvr^^^ = llgjM, , 

o dMo 9 

16^'^ = 52 M2 , 

16^^^ = -3ffiM3 (3) 



A. 3 Squark and slepton masses 



Let 



? = Tr[m^ - 2m^ + m|, - m| + m|] + mj^^ - m\j^ , 
= ^Q, + + m\j^ + A] , 
= ^Qs + ml, + mjj^ + Al , 
m2 = ml, + ml, + ml^ + Al , (4) 



where the terms in bold are matrices with respect to generations. Then the RG equations 
are then 
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16^2 ^f" = 25a3yr^' - 45?Mf + gU (5) 
A. 4 MSSM Higgs masses 



16vr2 — ^ = 3yiM^ + y^M^ - g^Mf _ Sgjul - -g^ (6) 



A. 5 MSSM Trilinear couplings (Rescaled) 



levr^^ = ^y^rAr + Syfe^Ab + 3gjMi + 3^1 M2 

16vr' -rr = ^VbAb + y'^r + SgfMi + 3g^M2 (7) 



~dt 

A. 6 The and Bfi terms 



levr^^ = ^(3y2 + 3y2 + y2 _ ^2 _ 3^2 > ^ 



16^2 ^ = (3y2 At + 3^,2^5 + A, - g^Mi - 3yiM2 ) . (8) 



Appendix B. 1-Loop NMSSM renormalisation group equations 

Again the RGEs [45, 46] are written assuming the running is dominated by the third family. 
We as before define t = log{Q^/Ql) and g^ = |(c?fuT)2_ r^^^ NMSSM superpotential and 
soft terms are as defined in Eqs. (1.6) and (1.7). 
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B.l Gauge and Yukawa couplings 



dt 


= ngf 


dt 


4 

= 92 , 


dt 





(9) 

le^r^ ^ = + + - 3^? - 3gi) (10) 



167r^ ^ = [Sy^ + Syf + + 4A^ + - 5? - 3^i) 



dK- 
~dt 

B.2 Gaugino masses 



16;r2^ = Ac2(6A2 + 6«2) (11) 



16.2^ = glM, , 

leTT^^ = (12) 



B.3 Squark and slepton masses 

Let 



e = Tr [ml - 2nil + mf, - mi + m|] + m]j^ - mjj^ , (13) 



= ^Q, + ml, + mjj^ + Al , 
Mr = ml, + ml, + ml^ + , 

Then the RG equations are 
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IGtt^^ = 26azylMl - ^-g\Ml - ^g^Ml + ^gf^ 
IStt^^ = SawlM^r - glMl - Sg^Ml - -g^ 
167r2^^ = 2Sa3ylM^ - Agfuf + g^ 

B.4 NMSSM Higgs masses 

Let 

Ml = ml^ + mjj^ +m% + Al, 
Ml = 3m% + Al , 



167r^^ = ^VbM^ + VrM^ + A'M| - gjMf - Sg^Ml - ^gf^ 

dt X K 

B.5 NMSSM Trilinear couplings (Rescaled) 

167r2^ = QytAt + y^Ab + X^Ax + + 3^1 M2 + y^l^^s 

IBtt'^ = 6^6^ A + y^At + y^Ar + X'Ax + ^gfM^ + 3gjM2 + y^l^a 

167r2^ = 4y2^^ + Sy^A^ + X^Ax + SgfM^ + S^^Ms 
dA 

I^tt'-^ = ^ylAh + y^A^ + x'^Ax + SgjMi + S^fMs 



= AX^Ax + 3y|^t + 3ylAb + y^A^ + 2^2^^ + gjMi + S^IMs 
= 6k^A^ + 6X^Ax 



B.6 Other parameters of the NMSSM 




(18) 



drun 



+ 2A^(3mi + 2^xAx) + 2XKm% - ff?(m| - 2/iMi) - 85! (mi - 2//M2) 



+ 2X{ml{Ax + fi') + n{mjj^ + m^J) + K(m'^(^« + fi') + 2/xm^) 

(20) 
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